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Abstract 

In [Hj, the author derived an exact rotation-strain model in two dimensions 
for the motion of incompressible viscoelastic materials via the polar decomposi- 
tion of the deformation tensor. Based on the rotation-strain model, the author 
constructed a family of large global classical solutions for the 2D incompressible 
viscoelasticity. To get such a global well-posedness result, the equation for the 
rotation angle was essential to explore the underlying weak dissipative structure 
of the whole viscoelastic system even though the momentum equation for the ve- 
locity field and the transport equation for the strain tensor have already formed 
a closed subsystem. In this paper, we revisit such a result without making use 
of the equation of the rotation angle. The proof relies on a new identity satisfied 
by the strain matrix. The smallness assumptions are only imposed on the H 2 
norm of initial velocity field and the initial strain matrix, which implies that 
the deformation tensor is allowed being away from the equilibrium of 2 in the 
maximum norm. 

Keyword: Rotation-strain model, viscoelastic fluids, partial dissipation, weak dissi- 
pative mechanism. 

1 Introduction 

In this paper, we revisit the two dimensional rotation-strain model for the incompress- 
ible viscoelastic fluid flows: 

V ■ u = 0, 

ut + u ■ Wu + = (iAu + V ■ (VV T ) + 2V ■ V, 

V t + u-VV = S(u) + \{WuV + VVu T ) (1.1) 
+ \u 12 {u){VA - AV) - \i{VA - AV), 

e t + u- w = -w 12 (u) + 7. 
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Here u(t, •) : M? — > M. 2 is the velocity field, p(t, •) : 1R 2 — > R is the scalar pressure, V 
and 6 are the strain tensor and the rotation angle of the deformation tensor, \i > is 
the viscosity coefficient and A is a constant matrix 

A=r?-A (i.2) 



x i o 7 

Besides, denotes the vorticity tensor: 

u ( u ) = ~ Vw T ), u 12 (u) = ^{d 2 ui - d x u 2 ) 

S(u) denotes the symmetric part of the gradient of velocity 

1 _ — T\ 



and 7 is given by 



S( u ) = -(Vu + Vu 1 ) 



7 = 2 + trV" ^ tVu3i2 ( u } ~ i^kUiV k2 - V k u 2 V k i)}. (1.3) 



System (II. ip was derived in [T3] from one of the most basic macroscopic models for 
viscoelastic flows 

\ -u = 0, 

u t + u ■ V"u + Vp = fiAu + V • (FF T ) , (i.4) 
^ + u • VF = VmF, 

via the rotation-strain decomposition of the deformation tensor 

F — (I + V)R, (1.5) 

where i? is orthogonal 

FF T = / 

and / + V is positive definite and symmetric 

V = V T . 



Physically, this means that the deformation is decomposed into stretching and rotation. 
Following a suggestion by K. O. Friedrichs [6], the tensor / + V is called the left 
stretch tensor, V the strain matrix and R the rotation matrix. We refer to [7J [23] for 
detailed discussions on the relevant physical background of system (11.41) and [H] for 
the derivation of (11. ip . 
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Note that the momentum equation for the velocity field u and the transport equa- 
tion for the strain tensor V in (11 .ip have already formed a closed subsystem: 

r u t + u ■ Vu + Vp = Am + V • {VV T ) + 2V • V, 

< V t + u-VV = S{u) + \{VuV + VVu T ) (L6) 
+ \u X2 {u){VA-AV)-\i{VA-AV), 
V -u = 0. 

We have set the viscosity fi to be one in the above system since we are not concerned 
the limit \i — > 0+ here. For the initial data 

u{0,x) =u {x), V{0,x) = V o {x), xeR 2 , (1.7) 

we impose the following constraints 

V ■ u = 0, 

det(/ + V Q ) = l, (1.8) 
V-V = A{I + V )V8 . 

The main result of this paper is the following theorem: 

Theorem 1.1. There exist two positive absolute constants C > 1 and e < 1 such 
that there exists a unique global classical solution (u, V) to the Cauchy problem for the 
rotation- strain viscoelastic model (jl.6p and (11.71) with the intrinsic physical constraints 
(11. 8p on the initial data provided that uq, Vq e H 2 and 

Ikoll^ + HVbllna < e o- 
Moreover, the solution (u, V) satisfies the following a priori estimate: 

POO 

^p{\Ht,.)\\ 2 H2 + \\V(t,-)\\ 2 H2 )+ / (HVulllp 

t>0 Jo 

+ \\VV\\ 2 Hl )dt < C(\\u \\ 2 H2 + \\V \\ 2 H i). 

Remark 1.2. 1). The constraints on the initial data in (11.81) are all the consequences 
of the incompressibility. See [2] for more details. 2). The result is also true in the two- 
dimensional torus case and bounded domain case. The periodic domain case can be 
treated similarly as in this paper. The bounded domain case is more involved. One can 
follow, for instance, the method in [21J. 3). Even though we don't utilize the equation 
for the rotation angle 8, there is still a smallness constraint on the norm ||V#o||l 2 - 
However, the smallness of ||V#o||l 2 does not imply that 8 is a small perturbation from 
a constant angle, and thus the deformation tensor F can be away from the identity at 
the distance of 2 in L°°: 

\\F(t, ■) - /|U~ = ||/ + V(t, •) - R T \\ L °° < 2. 
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This constraint on 9 (see the last equation in fll.il) ) is due to the inherent incompressible 
constraints (II. 8p . Besides, the dynamics of the rotation angle 6 is still inherently 
determined by the dynamics of velocity field and the strain tensor. 4). It is expected 
to study this problem in critical Besov space, but we do not pursue such a result in 
this paper. 5). Theorem 11.11 can be extended to the case of general energy function 
W = W(F) by the same argument in this paper. 

We emphasis that similar result in Theorem 11.11 has been proved in [H] by exploring 
the weak dissipative mechanics of the whole system (II. ip . where the proof essentially 
relies on the use of the equation for angle. We avoid using the equation for angle in 
this paper by deriving the following new structural identity (see Lemma 12. ip 

V • V • V = -V • [AV(I + V)- 1 AW ■ V] . 

This new identity together with the following Hodge decomposition (see step 4 of 
section 2 for its derivation) 

AV = VV • V + (V^trV - V ± (AV ■ V) 

is sufficient to prove that both u and it + 2// -1 A~ 1 V • V are dissipative using the similar 
arguments as in |14j . 

Let us end this introduction by mentioning some related works on incompressible 
viscoelastic system of Oldroyd-B type, which has attracted numerous attentions in 
the past and recent years. The study of the different contributions of strain and 
rotation can be traced back to Friedrich j6] where he observed that the smallness of 
the strain in nonlinear elasticity can be realized through the polar decomposition of 
the deformation tensor. John [TTj [T2] showed that no pointwise estimate for rotations 
in terms of strains can exist even in the case of small strain. In the work of Friedrich 
and John, no PDE is involved. From the PDE point of view, Liu and Walkington [23] 
considered some approximating systems resulting from the special linearization of the 
original system with respect to the strain. Lei, Liu and Zhou [16] constructed a 2D 
rotation-strain viscoelastic model and proved the global existence of classical solutions 
with small strain, but the dynamics of strain and rotation are not equivalent to that 
of the deformation tensor. The exact rotation-strain model was then derived in [16]. 

The global existence of classical solutions to ( 11 .4]) near equilibrium was established 
by Lin, Liu and Zhang [20J in the two-dimensional case (see also Lei and Zhou [19] 
via the incompressible limit method). The three dimensional case was then proved 
by Lei, Liu and Zhou [T7] (see also Chen and Zhang [3]). Very recently, an improved 
result was obtained by Kessenich [13J by removing the dependence of the smallness of 
the initial data on the viscosity via the hyperbolic energy method. We also mention 
the works on elastodynamics [U [29j [301 EB E2] and on viscoelastic Oldroyd-B model 

The paper is simply organized as follows. In section 2, after some preparations, we 
will present the proof of Theorem, which will be divided into four steps. 
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2 Proof of the Theorem 



Before the proof of Theorem II. 1[ let us give a brief explanation of the constraints 
on the initial data in (jl.8p . Due to the incompressibility, one has det F = 1. Since 
the deformation tensor F is decomposed into stretching and rotation F — (I + V)R, 
one has det(7 + V) — det F = 1, which is the second equality in (jl.8p . The third 
equality in (11. 8p can be derived from the identity V • F T = together with the polar 
decomposition of the deformation tensor. All of these constraints in (I1.8P are intrinsic 
for the viscoelastic Oldroyd model ( 11. 4ft . In fact, these identities in ( 11. 8ft are preserved 
in time (see [T4"]): 

fdet(I + V) = l, 

|v-y = A(/ + y)v6>. 1 ' J 

A direct consequence of the above intrinsic properties is that the linear terms and 
nonlinear terms appearing in the equations are not transparent any more, which is 
revealed in the following lemma: 

Lemma 2.1. Let F = (I + V)R be the polar decomposition and det F — 1. Then there 
hold 

trV = - det V (2.2) 

and 

V -V -V = -V ■ [AV(I + V)- 1 AV -V], (2.3) 
where A is given in ( 11. 2ft . 

The proof of Lemma [2.11 is straightforward. In fact, ( 12. 2 p is equivalent to the first 
identity in ( 12. ip . To see this, by the second identity in ( 12. ip . one has that 

V9 = -(I + V)- 1 AV -V (2.4) 

and 

v • v = Ave + Awe. (2.5) 

Apply the divergence operator to ( 12. 5ft . one has 

V ■ V ■ V = V ■ (AW 9). 

Then ( 12. 3D follows by plugging ( 12. 4 p into the above equality. 

The two identities ( 12. 2 p and ( 12. 3 p make system ( 11.61) fully dissipative. It is not 
necessary to use the equation for angle to complete the proof any more once one has 
those identities which are inherent in the fluid motion. Even though the proof of 
Theorem 11.11 is similar to that in [T3] , but the treatment of the term involving pressure 
in ( I2.15P is different. Besides, we need the space-time estimate not only for VAC/, but 
also for AC/. Here U is an auxiliary function 

U = u + 2/j,' 1 AW ■ V. 
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For a self-contained presentation, we will still present the whole proofs, but we will 
omit some of the similar parts. The proofs will be divided into four steps. In the first 
step we will show the basic energy law (I2.9p . However, the basic energy law is not 
enough to give an estimate of the low frequency part of the whole strain matrix, which 
will be complemented by an L 2 estimate. In the second step we will do higher order 
energy estimate for u and V. We then estimate the space-time L 2 norms of AU and 
V AU in the third step. In the last step we apply Hodge's decomposition (I2.28P and 
the structural identity (12. 3p to close the argument. 

Step 1. Basic Energy Law and an L 2 Energy Estimate. 

Let us first take the L 2 inner product of it-equation with u and F-equation in ( 11. 61) 
with V, and then add up the resulting equations to get 



-T- / (\ U \ 2 +\ V \ 2 ) dx+\\Vu 

2 at J R 2 

= -/ M . v( N!±g +p)dx 

JR 2 1 
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- / <V,\{VuV + VVu T ) + \u l2 {u){VA-AV)>dx 
l R 2 2 2 



+ / < u, 2V • V > + < V, S(u) > dx 
Jr 2 

+ / < u, V • (VV T ) > dx - [ <V, \i{VA - AV) > dx. 
On the other hand, the V^-equation also gives that 

~ / 2tiV dx = - [ u- VtiV dx+ [ trS(u) dx 

2 dt J R 2 J R 2 J R 2 

+ / tr(VuV + -(VA-AV)(co 12 (u)-j))dx. 
Jr 2 2 

Adding up the above two equations and using the divergence-free property of u, we 
obtain that 

1 d f 

0^7/ {M 2 + \V\ 2 + 2trV) dx+\\Vu\\ 2 L 2 (2.6) 

2 at J R 2 

= [ < V,-(VuV + VVu T ) + -co 12 (u)(VA- AV) > dx 
Jr 2 2 2 



+ / < u, 2V ■ V > + < V, S(u) > dx + / tr(VuV) dx 
Jr 2 Jr 2 

+ I < u, V • (VV T ) > dx - [ <V, \^{VA - AV) > dx 
Jr 2 Jr 2 2 

+ / tri((yA-AV)(wi2(tt)-7)) dx. 

JlB>2 Z 
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Let us compute the terms on the right hand side of (12.61) . First of all, by noting 
that V is symmetric, it is rather easy to see that 

/ < u, 2V • V > + < V, S(u) > dx + / tr(VuV) dx = 0. (2.7) 
Jr 2 Jr 2 

On the other hand, using the definitions in (11.21) and fll.31) . we have 

I < V, \{VuV + VVu T ) + \u l2 {u)(VA- AV) > dx (2.8) 
Jr 2 2 2 

+ / < u, V ■ (VV T ) > dx - I <V, -~f{VA - AV) > dx 
Jr 2 Jr. 2 2 

+ / trl-((VA-AV)(u 12 (u)- 7 ))dx 
Jr 2 1 

= [ <V, VuV > + < u, V ■ (VV T ) > dx 
Jr 2 

+ I -(u 12 (u) - 7) (tr(VA - AV)+ < V, V A - AV > ) dx 
Jr 2 2 

= 0. 

Then, combining ( 12. 6p . ( 12. 7p with (12. 8p . we have the basic energy law 

\i I (M 2 + \ V \ 2 + 2ilV ) dx + II Vw lli 2 = °- ( 2 - 9 ) 
2 ctr J R 2 

However, noting ( 12. 3p . one has 

\V\ 2 + 2trV = \V U - V 22 \ 2 + \V 12 + V 2l \ 2 . 

Thus, the basic energy law does not give an L 2 estimate of V. On the meantime, we 
have 

f°° 1 

J o l|v«||i 2 ^<-(|H| 2 L2 + ||y|| 2 i2 )- (2.io) 

This estimate will be crucial to derive the space-time L 2 estimate of AU and VAC/ in 
step 3. 

Next, let us derive the L 2 energy estimate. By taking the L 2 inner product of the 
first and second equations in (11.61) with u and 2V, and then adding up the resulting 
equations and using integration by parts, we obtain 

~ / (|^ + 2|Vf)d*+||V M ||! 2 (2.H) 
l at J R 2 

< \\Vu\\ L 2\\V\\ 2 Li 

< \\V\\ L 2\\VU\\ L 2\\VV\\ L 2, 



provided that 

lltr V^IIloo < 1, 
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which is guaranteed by (12. 2p and the smallness of ||V||#2: 

||tr V|U«- = ||detV|U~ < IIVHI^. 

Step 2. Higher Order Energy Estimates. 

In order to exploit the higher order energy estimates, we apply A to equations for 
u and V of system (II. II) . and then take L 2 inner product of the resulting equations 
with Au and 2AV, respectively, to yield 



1 d 
2dt 



|A< + 2\AV\ 2 dx + (Am, VAp) 

|2 



r \Av\ 2 

+ / u-V( l —^ + \AV\ 2 )dx 

= (AV, A [VuV + V Vu T + W12 (u) (VA - AV) 

- ^iVA - AV)] \ + (Au, AV ■ {VV T )) + (Au, 2 AV ■ V) 
+ (2AV, AS(u)) + (Au, AAu) - (Au, [A(u ■ Vw) - u ■ VAu] 

- (^2 AV, [A(u ■ W) - u ■ VAV] \ . 

By careful calculations, we have (for the details, see pH] ) 

~ f (|Aii| 2 + 2|AV| 2 )^ + ||VAn||| 2 ) (2.12) 
2 dt J R 2 

< (\\u\\ H 2 + \\V\\ H 2)(\\Vu\\ 2 L2 + \\VAu\\ 2 L2 + \\AV\\ 2 L z). 

Step 3. Weak Dissipation for the Strain Matrix V . 

Let us recall the definition of the auxiliary function U : 

U = u + 2A _1 V • V (2.13) 

Below we will estimate the space-time L 2 norm of both AU and VA£7. The part for 
f || V AU\\ 2 L2 ds has been done in [H]. However, to get the fully dissipative nature 
of the whole system (II. ip . the argument in [T3] involves another auxiliary function 
6 = Am + 2V" L 6 ) . We are going to avoid using 9 below, and thus provide the global 
well-posedness of solutions to the subsystem (II. 6p . 
The equation for U is as follows: 

AU t + u ■ VAU + VAp = A 2 U + Au + /, (2.14) 



S 



where 



/ = - [A(u -Vu)-u- VAu] - [V • (u ■ W) 

- u ■ V(V ■ V)] + AV ■ (VV T ) + V ■ [(VuV + VVu T ) 
+ u 12 (u)(VA - AV) - ^iVA - AV)] . 

Take the inner product of (I2.14p with AU in L 2 (IR 2 ), and use the similar arguments as 
in ( 12. 71) . we have 

~ / \AU\ 2 dx + \\VAU\\ 2 L2 = (Ap,V ■ AU) + (Au,AU) + (f,AU). (2.15) 
i at j R 2 



In [14] , the following estimate has been obtained: 

l(/,AC/)| < (\\u\\ H , + \\V\\ H ,)(\\Vuf H2 + ||VAC/||| 2 + \\AV\\l 2 ), (2.16) 

\(Au,AU)\<±\\VAUf L2 + \\Vu\\ 2 L2 . (2.17) 

Here we treat the pressure term by making use of the structural identity ( 12.31) in 
Lemma 12.11 We first apply the divergence operator to the momentum equation of 
system (11.61) to get 

Ap = -tr(VuVu) + V • [V • (VV T )} + 2V • (V • V). 

Plugging ( 12.31) into the above inequality, one has 

Ap = -V • V • (u g> u) + V • [V • (VV T )} - 2V ■ [AV(I + V)^AV ■ V] . (2.18) 

Consequently using (12.181) . we have 

|(Ap, V ■ AU)\ < (\\u\\ H 2 + \\V\\ H 2)(\\Au\\ 2 L2 + \\AV\\ 2 L2 + \\VAU\\ 2 L2 ). (2.19) 

The combination of (gZESJ , flZHSD , d2~T7j) andflOl gives 

|||A^||i 2 + ||VAf/||i 2 <||V M ||i 2 (2.20) 
+ (\\u\\ H 2 + ||F|| H2 )(||V«||^ + \\VAU\\ 2 L2 + ||AF||| 2 ). 

Next, let us derive a space-time L 2 estimate for AU. For this purpose, we first take 
the L 2 inner product of it-equation in (jl.6p with At/ to get 

\\AU\\ 2 L2 - [ AU-u t dx = I AU- [u- Vu + Vp- V- {VV T )]dx. (2.21) 

JR2 JIR2 



First of all, one has 



f AU ■ u t dx = ~\\Vu\\ 2 L2 - 2 [ (V-V)-u t dx. (2.22) 



Secondly, by using (I2.3p . it is easy to see that 



AU ■ Vpdx = -2 / pV ■ V • Vdx 

Ir 2 



(2.23) 



= 2 1 pV • [AV(I + V) -1 AV • V] dx 
Jr 2 

= -2 / Vp ■ [i4V(J + V) -1 AV • V] dx. 
Jr 2 

Consequently by (12.231) and (I2.18p . we have 

/ AU ■ [u ■ Vu + Vp - V • (VV T )]dx (2.24) 
Jr 2 

< \\AU\\#\\u ■ Vu - V • (^ T )|| L2 + || Vp|| i2 ||AV(/ + V)~ 1 AV ■ V\\ L2 



<-\\AUf L2 + \\u-Vu-V-{VV 



TM|2 

II 2 



+ |Uy(/ + y)"Mv-y 



|2 
lL 2 



<d|At/||£ 2 



(2.25) 



(2.26) 



\\u\\h 2 + \\V\M 2 {\\Vu\\1 2 + \\VV\\1 2 ). 
Now plugging (jZggp and fl2T2D into fl2T2B . we get 

||AC/||| 2 +|||V M ||i 2 -4^ 2 (V-y)-^ 

<(lhll^ + ll^b 2 ) 2 (liv^||i 2 + ||vy||i 2 )- 

Taking the L 2 inner product of V^-equation in (11.61) with 4Vw, one has 

-if uV-V t dx = 4 I Vu[-u- VV + S(u) + -(VuV + VVu T 
Jr 2 Jr 2 2 

+ ^co 12 (u)(VA - AV) - -"f(VA - AV)]dx 

< livtiiii, + (little + ||viM 2 (||v w ||i 2 + nwiii,). 

Adding up (12.251) and (12.261) . one gets the following: 

^||Vu||£ 2 + Hl! 2 
< (\Mh 2 + \\v\\ H *) 2 (\\Vu\\l 2 + ||w|| 2 i2 ) + HVtilli,. 

Step 4- Weakly Dissipative Mechanics of (II. 6p . 

This step ties everything together. Firstly, by Hodge's decomposition, we derive 
from straightforward calculations that 



(2.27) 



Al/ = VVT-VxVxl/ 

V 2 (ViVi2-V 2 Vii) -Vi(ViVi2-V 2 Vii; 
V 2 (ViF 22 - V 2 F 2 i) -Vi(ViF 22 - V 2 V 2 i] 

W • V + (V^hrV - V ± (AV ■ V). 



(2.28) 



VV -V 
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Thus, by Lemma [2. 1} we obtain 

IIWIIhi < ||V- VUffi + HVdetVllia 

< \\V-V\\ L 2 + \\V\\ H 2\\AV\\ L 2, 

which gives that 

HVVlffi < ||V-V|| H i, (2.29) 

provided that ||V||h2 is sufficiently small. Using ( 12. 13ft and ( 12.29ft . we have 

||VV||hi < \\&u\\ H i + \\AU\\ H i. (2.30) 

Plug ([CTD into (l2~TTj) . (12~T21) . (12T2DD and (l2~27j) . and then combine them to derive 
that 

4 / (h\\ 2 H 2 + \\V\\ 2 H 2 + \\AU\\l 2 )dx + (||V W || 2 H2 + || AU\\ 2 m ) (2.31) 
< (\M\e» + l|y||^)(||V«||^ + ||AC/||^) + ||V M || 2 i2 . 
Noting (I2.10p . one can easily conclude from (12.311) that the following a priori estimate 

su P {\\u(t,.)\\ 2 H2 + \\V(t,.)\\ 2 H2 ) + 

t>0 

+ ||(a m + 2v • v)\\ 2 H1 ]dt < c(\\u \\ 2 H2 + ||y ||| 2 )- 

is true if the H 2 norms of uq and Vq are sufficiently small. This further gives that 

)\(Au + 2V-V)\\ 2 m 



\W 2 



ff 2 



oo 



2 / \\V ■ V\\ 2 m dt < 



o Jo 

dt<C(\\u \\ 2 H2 + \\V \\ 



\Au\\ 2 m 



Then the a priori estimate in Theorem 11.11 follows by the above inequality and (12.291) . 
The proof of Theorem 11.11 is completed. 
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